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A NEW  METHOD  FOR  DESIGNING  SHOCK-FREE  TRANSONIC  CONFIGURATIONS 

t +■*■ 

H.  Sobieczky*,  N.  J.  Yu**,  K-Y.  Fung  , and  A.  R.  Seebass 

University  of  Arizona 
Tucson,  Arizona  85721 

Abstract 

A new  method  for  the  design  of  shock-free  supercritical  airfoils, 
wings,  and  three-dimensional  configurations  is  described.  Results  illus- 
trating this  procedure  in  two  and  three  dimensions  are  given.  They 
include  modifications  to  part  of  the  upper  surface  of  an  NACA  64A410  air- 
foil that  will  maintain  shock-free  flow  over  a range  of  Mach  numbers  for 
a fixed  lift  coefficient,  and  the  modifications  required  on  part  of  the 
upper  surface  of  a swept  wing  with  an  NACA  64A410  root  section  to  achieve 
shock-free  flow.  While  the  results  are  given  for  lnviscid  flow,  the  same 
procedures  can  be  employed  iteratively  with  a boundary  layer  calculation 
in  order  to  achieve  shock-free  viscous  designs.  With  a shock-free  pressure 
field  the  boundary  layer  calculation  will  be  reliable  and  not  complicated 
by  the  difficulties  of  shock-wave  boundary-layer  interaction. 
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INTRODUCTION 


Well-known  requirements  for  increased  efficiency,  and  in  the  case  of 
aircraft,  productivity,  have  forced  the  operating  conditions  of  compres- 
sors, turbines,  propellers,  wing  sections,  and  aircraft  into  the  transonic 
regime.  Unfortunately,  once  local  regions  of  supersonic  flow  occur,  shock 
waves  are  likely  with  the  attendant  wave  drag,  and  boundary  layer  separa- 
tion, losses.  In  the  mid-fifties,  Morawetz1  proved  that  shock-free,  two- 
dimensional,  irro rational,  near  sonic  flows  are  mathematically  isolated. 

In  other  words,  any  small  changes  in  the  flow  or  boundary  conditions  that 
provide  a shock-free  flow  will  lead  to  the  formation  of  a shock  wave.  Thus 
Morawetz's  theorem  stated  that  the  shock-free  invlscid  flow  solutions,  if 
and  when  they  existed,  were  Isolated  by  neighboring  solutions  that  contain 

shock  waves.  Recently  this  result  has  been  extended  to  three  dimensions 
2 

by  Cook.  Fortunately,  it  was  recognized  that  such  flows  would  have  prac- 
tical significance  if,  as  seemed  likely,  the  shock  waves  that  occurred  in 

3 

neighboring  flows  were  very  weak.  Wind  tunnel  research  by  R.  T.  Whitcomb 

at  the  NASA  Langley  Research  Center  and  by  H.  H.  Pearcey  at  the  National 

Physical  Laboratory  (U.K. ) led  to  the  development  of  practical  "shock-free" 

airfoil  sections.  Subsequent  analytical  studies  by  Garabedian  and  Korn,^ 
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Nieuwland,  Boerstoel,  and  Sobleczky  established  theoretical  design 
procedures  for  two-dimensional  invlscid  flows.  More  recently,  the  develop- 
ment of  sophisticated  numerical  codes  for  the  analysis  of  transonic  flow 

fields  has  led  to  the  design  of  both  airfoils  and  wings  by  numerical 
9,10 

optimization.  The  practical  success  of  the  above  efforts,  as  documented 

by  the  recent  NASA  Conference  on  Advanced  Technology  Research, ^ has  been 
substantial.  Further  progress,  as  reported  here,  seems  likely.  The  senior 
author  recognized  that  the  procedure  he  was  using  in  the  hodograph  plane 


Implied  an  analogous  procedure  in  the  physical  plane,  and  further,  that 
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this  procedure  did  not  seem  to  be  restricted  to  two-dimensional  flows.  * 
This  paper  reports  the  success  we  have  had  to  date  in  using  this  idea  to 
provide  shock-free  designs  in  two  and  three  dimensions. 

The  design  procedure  invoked  here  is,  in  principle,  a simple  one. 

While  there  is  no  guarantee  that  a shock-free  flow  will  necessarily  result 
from  the  procedure,  our  experience  in  two -dimens ions  has  been  chat  if  the 
holograph  method  will  work  for  specified  flow  and  airfoil  parameters,  then 
the  procedure  outlined  here  will  work,  coo.  Also,  it  provides  neighboring 
shock-free  airfoil  shapes  for  fixed  lift  coefficient  with  varying  Mach 
numbers  and  varying  lift  coefficient  for  fixed  Mach  numbers,  as  well  as 
providing  a multiplicity  of  closely  related  shapes  that  are  shock-free  at 
fixed  lift  coefficient  and  Mach  number.  This  wealth  of  shock-free 
two-dimensional  designs  is  of  no  great  surprise;  It  is,  Chen,  not  surpris- 
ing that  they  are  found  with  minimal,  computational  effort.  Two-dimensional 
invisc id  flow  potential  airfoil  designs  require  less  chan  a minute  of 
CYBER  175  CPU  time  and  only  a few  seconds  of  CSC  7600  CPU  time. 

For  three-dimensional  flows  our  results  are  less  extensive.  Also, 
while  it  is  clear  chat  the  procedure  we  use  rests  on  a sound  mathematical 
foundation  in  two  dimensions,  this  may  not  be  the  case  in  three  dimensions. 
Indeed,  for  three-dimensional  (chat  Is  non- planar  and  non-axiayrmetrlc ) 
flows  we  are  probably  solving  an  Ill-posed  boundary  value  problem.  The 
fact  that  shock- free  flows  are  obtained  in  the  cases  studied  here  are  a 
consequence  of  the  pseudo-analytic  character  of  the  initial  data  and  the 
particular  numerical  technique  used  to  calculate  the  flow  in  the  hyperbolic 
region.* 


*The  authors  are  indebted  to  Professor  A.  Jameson  of  the  Courant  Institute 
for  alerting  them  to  this  difficulty. 
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We  have  demons traced  the  ability  to  modify  three-dimensional  wings 

so  that,  within  the  context  of  the  numerical  algorithm  used,  shock-free 

flows  are  obtained.  We  have  not  yet  demonstrated  an  analogous  wealth  of 

shock-free  flows  in  the  three-dimensional  case,  but  see  no  reason  to 

believe  that  this  situation  is  different  there.  The  practical  consequences 

14 

of  this  wealth  should  prove  to  be  of  interest  to  the  aircraft  industry. 

DESIGN  PROCEDURE 

The  procedure  we  use  to  find  shock-free  designs  assumes  that  a reliable 
numerical  code  is  available  for  computing  the  flow  past  a given  configuration, 
such  as  that  sketched  in  Figure  1.  Such  codes  are  available  for  two-  and 
three-dimensional  iavlscld  flows.  When  they  are  coupled  with  a reliable 
boundary  layer  code,  the  design  procedure  outlined  here  can  be  used  to 
calculate  shock-free  viscous  flow  designs.  While  this  would  require  some 
modest  iteration.  It  is  certainly  possible,  both  in  practice  and  in  prin- 
ciple. With  the  existence  of  a reliable  analysis  algorithm  presumed,  we 
modify  this  algorithm  so  that  once  the  flow  becomes  hyperbolic  we  alter  the 
basic  equations  so  chat  they  revert  to  elliptic  behavior.  This  may  be  done 


in  a number  of  ways,  but  it  should  be  done  in  a way  chat  it  conserves  new. 


but  fictitious,  "mass"  and  "momentum"  fluxes  to  a satisfactory  degree  of 
accuracy.  We  may,  for  example,  change  the  density's  dependence  from  the 
usual  one  to  one  that  returns  the  equations  to  elliptic  form.  We  might 
suppose,  for  the  purpose  of  illustration,  chat  once  the  aquations  become 
parabolic,  l.e.,  sonic,  on  some  surface  Chen  at  higher  velocities  the 
density  will  be  maintained  at  its  sonic  value,  giving  elliptic  equadons. 
We  use  a numerical  algorithm  to  compute  this  fictitious  flow  past  a con- 
figuration of  interest,  chosen  perhaps  on  the  basis  of  previous  design 


I 
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experience.  Because  the  aquations  are  elliptic  this  will  result  la  a 
discretized,  pseudo-analytic,  description  of  the  velocity,  density,  and 
pressure  fields  on  the  embedded  parabolic  surfaces,  and  this  description 

1 

will  be  consistent  with  the  correct  governing  equations.  This  initial 

data  on  the  parabolic  surfaces  is  then  used  to  calculate  the  correct  flow 

field  inside  such  surfaces.  This  new  flow  field  may,  or  may  not,  contain 

shock  waves.  This  depends  on  the  choice  of  the  ficticious  equations,  or 

perhaps  better,  fictitious  gas,  used  inside  the  parabolic  surfaces.  This 

Q<v  flov  will  define  a scream  surface  chat  is  tangent  to,  has  the  same 

t\ 

curvature  as,  the  stream  surface  at  the  Intersection  of  the  sonic  surface 
and  the  original  body.  Inside  this  surface  a new  body  shape  is  defined  by 
the  stream  surface  of  the  new,  but  now  real,  flow. 

Here,  of  course,  we  must  also  address  the  question  of  whether  or  not 
this  initial  value  problem  is  well  posed.  In  two  dimensions  there  is  no 
difficulty  because  either  of  the  spatial  coordinates  may  be  designated  as 
the  time— like  variable.  This  is  not  the  case  in  three  dimensions  where 
only  the  spatial  coordinate  aligned  with  the  flow  is  time-like.  Because 
shock-free  flows  are  reversible,  the  domains  of  dependence  aui  influence 
may  be  interchanged.  But  neither  the  normal  (nor  the  binormal)  to  the 
scream  direction  can  be  considered  time-like  in  the  three-dimensional 
initial  value  problem.  Thus,  it  may  be  ill-posed  because  data  are  given 

I 

on  surfaces  chat  are  not  in  the  uaual  domain  of  dependence.  If  so,  any 
computational  algorithm  will  be  unstable  for  the  three-dimensional  problem. 

I ! Further,  while  such  computations  can  be  stabilized  by  artificial  means,  the 

* 

results  must  be  considered  suspect  until  they  are  verified  by  an  independent 
computation.  It  is  this  fact  chat  has  made  us  stress  that  a reliable 

algorithm  should  be  the  basis  for  the  design  computations.  For 

■I 
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tvo-dimens local  (planar  and  axisymmetrlc)  designs  this  difficulty  does  not 
occur  because  the  lateral  coordinate  can  be  considered  to  be  the  time-like 
direction.  A simplistic  analysis  of  model  problems  Indicates  that  varia- 
tlons  In  the  spanwise  direction  that  are  on  a scale  that  Is  small  compared 
to  the  nominal  axial  (flow  direction)  distance  may  amplify;  thus  the 
success  of  the  numerical  algorithm  here  may  depend  upon  Its  natural  filter- 
ing of  such  disturbances.  This  Is  not  the  first  time  Ul-posed  problems 
have  been  solved  to  obtain  results  of  engineering  Interest;  see,  for 
example.  Ref.  15,  pp.  448-472. 

Fictitious  Gas 

As  mentioned  above,  modifications  are  made  to  the  basic  equations  to 
retain  their  elliptic  behavior  once  the  flow  has  accelerated  to  sonic  speed 
and  a parabolic  surface,  with  the  needed  initial  data,  has  been  generated. 

The  possible  modifications  are  manifold.  We  limit  our  discussion  to  those 
we  have  used  to  obtain  the  results  reported  here. 

For  two-dimensional  flows  we  have  used  Jameson ' s^*^  circle-plane 
algorithm  for  the  full  potential  equation.  Thus,  in  the  analysis  mode,  we 
are  solving 

+ (p$  } - 0 (la) 

• 4 z z 

with 

p/pm  - (1  + 1-y-ii£(l  - *1  ~ ^)]1/Y-1  (lb) 


where  <p  is  the  velocity  potential  and  p the  density.  If  we  limit  our 
consideration  to  fictitious  gases  for  which  the  density  is  a function  of 

l 


2 2 2 2 2 

the  square  of  the  velocity,  viz.,  p - p (q  ),  where  q - 0 [$“  + $al, 
then  gas  laws  of  the  form 

p 

P /p*  - U*/q)  » P < 1,  for  q > a*  (lc) 


will  insure  elliptic  behavior;  P • 1 gives  parabolic  behavior  and  the 
fictitious  and  real  gases  have  the  same  value  of  (dp/dq)A.  An  alternative 
choice,  and  the  one  we  have  used  most  extensively  here,  is  P ■ 0;  in 
this  case  Equation  (la)  becomes  Laplace's  equation.  When  the  flow  would 
normally  be  hyperbolic  we  now  solve  Equation  (la)  with  the  density- 
velocity  relationship  of  Equation  (lc) . A fictitious  mass  flow,  which 
matches  the  real  mass  flux  at  the  sonic  surface,  is  thereby  conserved  and 
the  velocity  field  remains  irro rational. 

For  three-dimensional  flows  we  have  used  the  Ballhaus,  Bailey,  Frick 
18  19 

algorithm,  as  implemented  by  Mason  et  al.  This  is  a small  perturba- 
tion calculation  and  we  adopt  the  classical  conservative  formulation  here. 
Thus  we  solve,  in  an  equivalent  form,  the  system 


- -r  (y  + lHu2}  +>  v + w - 0 
i.  x y z 


u 


7 


v 


X 


0 


(2) 


0, 


where  the  velocity  vector  is  <j^  ■ a*[(l  +■  uH  + vj_  + wk] . 

2 2 

A simple  modification  (2)  is  to  replace  (u  } by  -sgn(u){u  } 

X X 

for  all  u.  This  system  is  elliptic  except  on  the  sonic  surface  where 


u - We  may  chink  of  the  first  of  Equations  (2)  as  being  the  consequence 
of  the  snail  perturbation  expansion  for  the  density,  viz.. 


(J  . _ y - 1 2 

_ 1 ■ -u  - = — u . 


whereas  the  fictitious  equation,  with  u replaced  by  -|u|  for  u > 0, 
results  from 


f-  - 1 - -»  + X-t-l  u2; 

0*  2 


this  fictitious  gas  has  the  same  value  for  (do/duK  as  the  real  gas. 


Equation  (3).  For  three-dimensional  design  studies,  then,  we  solve  Equa— 


2 2 

dons  (2)  with  {u  replaced  by  -sgn(u){u  this  corresponds  to 


using  the  densities  given  by  Equations  (3)  and  (4)  for  u < 0 and  u > 0 
respectively. 


Calculation  of  the  Hyperbolic  Flow  Field 


As  described  above,  we  calculate  the  flow  past  a body  using  the  correct 
equations  when  the  flow  is  subsonic  and  a modified,  incorrect,  set  of  equa- 
tions when  the  flow  is  supersonic.  This  calculation  serves  to  define  sonic 
surfaces  on  which  the  flow  field  calculadon  is  switched  from  the  correct 
•^wations  to  the  modified  ones.  Outside  this  surface,  presuming  the 
trailing  edge  of  the  wing  is  subsonic,  the  solution  satisfies  the  correct 
equations  and  the  potential  at  infinity  has  the  correct  value  for  the 
circulation.  If  infinity  in  the  physical  plane  is  not  mapped  to  a finite 
P*rc  che  computational  plane,  then  there  is,  in  principle,  a need  to 
correct  the  doublet  and  nonlinear  contributions;  in  practice,  these  con— 


r 

cribu  cions  are  small  and  changes  in  chon  negligi-xe.  Thus  che  flow  in  the 
elliptic,  subsonic,  domains  is  fixed  and  known,  as  is  che  initial  data  we 
need  on  che  parabolic  surface. 

For  two-dimensional  flows  che  calculation  of  the  correct  hyperbolic 
behavior  is  carried  out  using  the  method  of  characteristics.  This  is  done 
in  a hodograph-like  working  plane  in  which  che  characteristics  are 
orthogonal  straight  lines.  If  we  take  £ ■ 9 + v and  n ■ 9 — \»  where 
9 is  the  flow  deflection  angle  and  \j  the  Prandtl-Meyer  turning  angle, 
Chen  che  velocity  potential  and  scream  function  satisfy 


h - 


(5) 


♦n  * 


or  equivalently. 


d1» 

d* 


- ±K 


-1 


C,n-consc 


where  the  t signs  refer  to  • const.,  respectively.  Here 


K(v)  - K[v(q)}  - { | M2(q)  - l|  }1/20  (0)/p  (q) 


Values  for  che  velocity  potential  on  the  parabolic  line,  z - z*(x),  and 
the  shape  of  this  line  are  used  along  with  che  usual  relations  between  the 
spatial  coordinates  and  $ and  1>  to  find  1»  on  the  sonic  line.  This 
initial  data  is  then  integrated  u^-og  Equations  (S)  to  determine  che  locus 
l/(x,z)  - 0 which  passes  through  cae  intersection  of  che  sonic  line  with 
che  body  surface.  The  values  of  z for  which  4i(x,z)  - 0 determine  the 


J 
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new  body  shape.  This  shape  will  have  che  same  slope,  and  at  lease  theo- 
retically, the  same  curvature,  as  the  original  body  at  the  sonic  points. 
This  follows  from  the  observation  that  flow  quantities  are  not  changed  at 
the  sonic  line;  thus  Che  streamwise  momentum  and  normal  pressure  gradient 
are  unchanged.  Consequently  the  local  flow  curvature  must  be  the  same. 

For  three-dimensional  flows  the  calculation  of  the  hyperbolic  flow 
field  is  carried  out  by  a procedure  that  marches  inward  from  the  sonic 
surface  by  successive  surfaces  of  constant  density  (isopycnics)  for  the 
full  potential  equation,  or  constant  axial  flow  speed,  u,  for  the  small 
perturbation  equation.  We  limit  our  discussion  to  the  small  perturbation 
equations,  as  all  the  results  reported  here  derive  from  them.  Preliminary 
results  using  the  full  potential  equation  have  been  obtained  by  one  of  the 
authors  (N.  J.  Yu). 

We  may  either  write  the  Equations  (2)  in  the  appropriately  scaled 
form  or  work  with  them  directly,  which  we  will  do  here. 

We  are  given  an  isotach  surface  r*(x,y),  as  shown  in  Figure  2,  on 
which  we  know  u • u*  • const.,  w - w*(x,y),  and  v - v*(x,y).  We  use 
the  data  on  this  surface,  and  the  surface  shape,  to  calculate 


z*,  z*,  w*  w* . v*,  v*. 

x y x’  y*  x’  y 


(6) 


Because  this  data  satisfies  Equations  (2)  we  can  verify  that 


v*  ■ z*w*  - 2*W* 
x x y y x 


which  can  be  used,  if  needed,  to  check  the  consistency  of  the  initial  data. 
The  values  given  in  Equation  (6)  can  now  be  used  to  calculate  the  z 


'! 


derivatives  of  u,  w,  v on  z*(x,y),  where  u(x,y,z*) 
using 


const.,  by 


u “ [z*v*  - z*v*  - w*l/J 
2 y x x y X 

w • [ (y  + l)u*z*v*  - z*v*  + v*]/J  (7) 

z x x y y yJ  ' ' 

v • [(y  +>  l)u*z*v*  - w*  - z*v*]/J 
2 x x y y y 

where  J,  the  Jacobian,  3 (u,v,w)/3 (x,y,z) , is 

J • (y  + l)u*z*2  - z*2  - 1. 
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When  the  Jacobian,  which  is  initially  negative,  vanishes  we  can  no  longer 
compute  the  z derivatives;  this  corresponds  to  the  subsequent  formation 
of  multi-valued  solutions,  i.e.,  limit  surfaces.  If  J ■ 0 occurs  before 
the  calculations  produce  a suitable  stream  surface  defined  by  w(x,y,0) , 
v(x,y,0),  then  they  must  be  rejected. 

With  the  first  of  Equations  (7)  inverted  to  give  (dz/du)A,  we  take  a 
sat  increment  in  u,  Au,  co  form  a new  isotach  surface  z*(x,y)  +■  lz*(x,y). 
This  new  shape,  along  with  the  mean  value  of  u between  the  two  surfaces 
and  the  second  and  third  of  Equations  (7),  provides  che  new  values, 
w*(x,y)  + iw*(x,y),  v*(x,y)  +■  Av*(x,y),  of  w*  and  v*  on  the  next 

isotach.  These  values  and  che  shape  of  che  subsequent  isotach  are  then 
converted  co  continuous  functions  by  one-dimensional  cubic  splines  in  the 
x and  y coordinates.  This  "onion-p«el,,-like  process  is  then  continued 
until  z - 0,  unless  a limit  surface  intervenes.  In  che  latter  tvent  che 
solution  must  be  rejected.  A more  detailed  discussion  of  chis  procedure 
is  given  in  Ref.  20. 
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TWO-DIMENSIONAL  RESULTS 

We  have  explored,  rather  extensively,  some  of  the  modif ications  that 
can  be  made  to  an  existing  airfoil,  namely  an  NACA  64A410  airfoil,  to 
obtain  shock-free  flow.  We  will  call  this  the  baseline  airfoil,  as  the 
airfoil  shapes  we  generate  are  identical  with  this  airfoil  over  that  por- 
tion wetted  by  subsonic  flow;  we  need  only  modify  the  airfoil  over  a 
limited  portion  of  its  upper  surface  to  obtain  shock-free  flows.  Further, 
this  modification  is  not  unique  for  fixed  flight  conditions;  rather,  if 
one  such  shape  exists,  there  will  be  an  infinite  family  of  modifications 
of  the  baseline  airfoil  that  will  produce  shock-free  flow. 

With  a baseline  airfoil  selected,  here  mainly  for  illustrative  pur- 
poses, we  then  pick  a sec  of  flight  conditions  for  which  we  wish  to  find  a 
modification  of  the  airfoil  shape  chat  will  result  in  shock-free  flow.  We 
choose  Ma  * 0.72  and  a,  the  angle  of  attack,  0.4  DEC.  At  these  condi- 
tions in  viscid  flow  calculations  for  the  NACA  64A410  baseline  airfoil  give 
a of  0.78  and  a of  0.0064.  The  design  procedure  discussed 

above  results  in  an  airfoil  that  is  9.3Z  thick  and  has  a lift  coefficient 
of  0.703.  The  original  and  the  design  pressure  coefficient,  sonic  lines, 
and  body  shapes  are  compared  in  Figure  3a;  these  results,  and  all  other 
"analysis"  results  were  computed  using  the  numerical  algorithm  of  Ref.  16. 
Figure  3b  compares  the  pressure  coefficients  and  sonic  lines  determined  by 
the  design  procedure  with  those  computed  for  the  design  airfoil  shape. 

With  this  shock-free  design  established  at  M ■ 0.72  and  with  CT  • 

la 

0.70,  we  now  wish  to  determine  the  families  of  shapes  that  provide  shock- 
free  flow  for  fixed  lift  coefficient  as  the  Mach  number  varies,  and  fixed 
Mach  number  a a the  lift  coefficient  varies.  This  we  have  done  with  P * 0, 
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chat  is,  with  a constant  density  fictitious  gas  (at  the  critical  value). 

We  have  then  explored  other  shapes  chat  vlll  produce  the  same  lift  coeffi- 
cient, 0.70,  at  a fixed  Mach  number,  for  three  different  Mach  numbers, 
by  caking  P to  be  -0.5,  0.5,  and  1.0.  Also,  for  P - 0 we  have 

determined  the  Mach  number  for  which  the  design  procedure  will 

produce  a shock-free  airfoil,  as  a function  of  lift  coefficient.  This 
Mach  number  is  nearly  a linear  function  of  Lift  coefficient  at  larger  lift 

coefficients.  The  slope  of  this  variation  is  consistent  with  chat  given  by 
21 

Boerscoel.  Preliminary  studies  also  indicate  chat  for  a fixed  lift 
coefficient  of  0.6-0. 7,  an  0.12  increase  in  Che  maximum  Mach  number 
requires  about  an  0.22  reduction  in  the  thickness  for  shock-free  flow, 
when  the  nominal  thickness  is  about  102.  This  result  is  less  optimistic 
than  the  envelope  of  the  hodograph  designs  given  by  Boerstoel, who  found 
chat  only  an  0.12  reduction  was  required.  In  our  study  the  generic 
family  of  the  airfoil  is  invariant;  we  have  not  yet  examined  the  modifica- 
tions required  when  the  baseline  airfoil  is  near  the  envelope  of  hodograph 
designs.  Positive  values  of  P provide  less  airfoil  thickness  reduction, 
as  the  ficticious  and  real  gas  densities  are  more  nearly  the  same.  The 
range  of  our  airfoil  studies  is  depicted  in  Figure  4,  with  shock-free 
airfoils  being  determined  for  the  points  indicated.  Also  shown  in  Figure 
4 is  the  maxinnm  Mach  number  for  which  a design  was  found  as  a function  of 
lift  coefficient  for  P ■ 0. 

The  accuracy  of  the  design  procedure  was  studied  at  a number  of  design 
points  by  comparing  the  design's  pressure  distribution  and  sonic  line  shape 
with  chose  obtained  using  the  unmodified  numerical  algorithm  to  analyze  the 
design  airfoil  shape.  Typical  results  are  shown  in  Figure  5.  The  sonic 
line  shape  and  initial  data  on  the  sonic  line  are  determined  in  the  circle- 


plane;  they  Chen  are  mapped  back  Co  Che  physical  plane.  The  mechod  of 

characcerlsdcs  in  the  hodograph  variables  Is  used  Co  compute  Che  design 

pressure  coeiflcienc  corresponding  Co  Che  calculated  airfoil  surface  shape. 

The  agreement,  as  shown.  Is  excellent.  For  designs  that  approach  the  Mach 

number  ac  which  a limit  line  firsc  penecraces  Che  surface  special  care  must 

be  taken  with  Che  analysis  code  In  order  co  obtain  a converged  soluclon. 

These  designs  have  very  rapid  expansions  immediacely  following  Che  sonic 

21 

line.  Indeed,  as  Boerscoel  has  no  ted,  che  analysis  code  used  with  an 
optimization  scheme  will  not  produce  designs  of  this  character. 

The  shock-free  airfoil  shapes  chat  are  obtained  for  fixed  and  P, 

fixed  M and  P,  and  fixed  M and  CT  at  various  P's,  are  shown  In 
Figures  6-6.  One  can  overlay  che  results  for  fixed  and  find  quite 

similar  airfoil  shapes  that  are  shock-free  over  a range  of  Mach  numbers. 
Because  modifications  to  che  baseline  airfoil  are  required  only  over  a 
limited  portion  of  che  upper  surface,  and  a family  of  specif led*  changes  In 
che  airfoil  curvature  Is  known  for  each  sec  of  flight  conditions,  a closely 
related  family  of  shock- free  airfoil  shapes  can  be  generated.  Thus  che 
minor  modifications,  co  a limited  portion  of  a wing  surface,  needed  co 
produce  shock-free  flow  over  a practical  range  of  flight  conditions  can 
easily  be  determined. 

THREE-DIMENSIONAL  RESULTS 

Our  first  design  results  using  the  mechod  described  above  were  for  two- 
dimensional,  small  perturbation  flow  past  a parabolic  arc  airfoil.  Conse- 
quently, we  Initiated  our  three-dimensional  studies  with  a rectangular, 
unswept  wing  with  an  aspect  ratio  of  six  and  a parabolic  arc  airfoil.  We 
utilized  che  small  perturbation  approximation.  Equations  (2),  and  a parabolic 


thickness  distribution;  the  airfoil  was  taken  to  be  6Z  thick  at  the  center 
plane.  The  flow  was  calculated  using  the  algorithm  of  Ref.  19,  modified  to 
return  the  equations  to  elliptic  behavior  as  described  earlier.  The  initial 
data  on  the  embedded  sonic  surface  was  then  used  to  compute  the  correct  flow 
in  the  supersonic  domain  using  the  "onion  peel"  algorithm  of  Ref.  20.  This 
defines  new  wing  surface  slopes.  The  flow  past  this  shock-free  design  was 
Chen  analyzed,  using  the  unmodified  numerical  algorithm.  Figure  9 compares 
the  pressure  distributions  on  the  original  and  design  wing,  at  various 
lateral  positions,  for  Ma  ■ 0.87.  Also  shown  are  the  cross  sections  of  the 
sonic  surface  at  the  same  lateral  stations.  The  only  essential  differences 
in  the  pressure  occur  in  the  supersonic  domain,  which  is  consistent  with  the 
design  process.  The  modifications  made  to  the  wing  slope,  shown  in  Figure 
10  for  several  lateral  stations,  have  eliminated  the  shock  wave. 

A subsequent,  more  realistic,  calculation  was  made  for  the  planform 
sketched  in  Figure  11.  The  wing  section  chosen  was  an  MACA  64A410  pro- 
file at  the  center  section  and  an  elliptic  thickness  distribution.  The 
leading  edge  sweep  was  taken  to  be  30  DEG,  the  trailing  edge  15  DEG  and 
the  span  to  chord  ratio  five.  The  sonic  surface  is  also  depicted  in  Figure 
11.  Figure  12  compares  the  pressure  coefficients  on  the  upper  surface  i£ 
the  original  wing  and  the  wing  designed  to  be  shock  free.  While  the  reduc- 
tion in  drag  for  this  wing  is  small  compared  to  the  induced  drag,  it  is 
clear  that  the  wing  modifications  have  essentially  eliminated  the  shock 
waves,  and,  consequently,  the  wave  drag.  More  importantly,  shock  wave 
induced  boundary  layer  separation  is  avoided. 

We  pause  at  this  point  to  stress  that  the  above  comparison  is  obtained 
by  computing  the  flow  past  the  original  wing  and  the  design  wing,  using  the 
same  numerical  algorithm.  The  process  that  leads  to  the  new  wing  shape 
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also  provides  Che  pressure  on  the  wing.  Because  chis  pressure,  and  the  new 
design  shape,  result  from  a problem  chac  is  presumably  ill  posed,  che 
results  were  considered  suspect  until  verified  by  che  original  numerical 
algorithm.  Our  experience  with  che  calculation  of  che  three-dimensional 
supersonic  flow  field  is  limited.  But  this  limited  experience  seems  to 
confirm  that  reducing  che  lsotach  or  lateral  mesh  size  to  small  values  gives 
results  that  are  indicative  of  an  Instability.  Ac  Chis  juncture  we  can  only 
say  that  results  of  engineering  interest  are  obtained  in  three-dimensional 
flows,  no  doubt  because  the  flows  of  interest  are  frequently  chose  whose 
behavior  is  generally  smooth  and  locally  analytic . 

CONCLUSION 

A novel  and  simple  procedure  for  determining  modif ications  that  will 
make  a baseline  configuration  shock-free  for  supercritical  flight  conditions 
has  been  delineated..  For  two-dimensional,  inviscid  flows,  shock-free  de- 
signs are  obtained  in  seconds  on  a CYBER.  175.  Families  of  airfoils  that  are 
shock-free  at  fixed,  as  well  as  varying,  flight  conditions  are  found.  The 
same  procedure  has  been  applied  to  three-dimensional  wings,  resulting  in 
wing  modifications  that  make  the  wings  shock-free  when  the  flow  is  analyzed 
with  che  numerical  algorithm  that  was  modified  to  become  a design  tool.  It 
can  also  be  applied  to  che  design  of  shock-free  cascades.  A unique  feature 
of  che  procedure  is  that  any  code  chac  is  effective  in  computing  che  flow 
field  may  be  modified,  in  various  ways,  to  be  a design  algorithm  if  it  is 
coupled  with  a method  for  calculating  the  solution  in  che  supersonic 
domains  for  given  data  on  che  sonic  surfaces.  A straightforward  marching 
technique  for  such  computations  is  described  for  three-dimensional  flows; 
in  two  dimensions  either  the  marching  procedure  or  the  method  of 
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charactariatics  may  be  used  for  tha  supersonic  domain.  The  algorithm  for 
tha  supersonic  domain  serves  to  define  the  modifications  needed  in  the 
conf iguration  to  achieve  shock-free  flow;  these  modifications  will  be 
limited  to  that  portion  of  the  design  shape  that  are  wetted  by  supersonic 
flow. 
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Skecch  of  shock-free  flow  pasc  a lifting  wing  depicting  the 
sonic  surface  obtained  by  introducing  ficticious  behavior 
inside  this  surface  chat  results  in  elliptic  equations.  The 
correct  flow  in  this  supersonic  domain  is  subsequently  calcu- 
lated using  the  initial  data  on  the  sonic  surface.  This 
calculation  provides  the  wing  geometry  modifications  needed 
to  obtain  shock-free  flow. 

Sketch  of  two  neighboring  iso each  surfaces  used  in  the  calcu- 
lation of  the  supersonic  domain  for  Equations  (2) . 

Comparison  of  the  pressure  coefficients  and  sonic  lines  for 
Che  baseline  NACA  64A410  and  the  shock-free  airfoil  obtained 
from  it  by  the  direct  design  procedure. 

Parameter  space  explored  for  the  shock-free  airfoils  that  can 
be  obtained  when  the  baseline  configuration  is  an  NACA  64A410 
airfoil. 

Comparison  of  the  pressure  coefficient  and  the  sonic  line 
obtained  by  the  design  calculation  that  modifies  the  airfoil 
shape  with  those  obtained  by  computing  the  flow  past  the 
modified  airfoil  with  the  numerical  algorithm  of  Ref.  16. 

Shock-free  airfoil  shapes  for  fixed  lift  coefficient  and  vary- 
ing Mach  number.  The  fictitious  gas  has  a constant  density  in 
the  supersonic  domain.  The  vertical  scale  is  magnified  five 
times  and  the  baseline  airfoil  is  an  NACA  64A410. 

Shock-free  airfoil  shapes  for  fixed  Mach  number  and  varying 
lift  coefficient.  The  fictitious  gas  has  a constant  density  in 
the  supersonic  domain.  The  vertical  scale  is  magnified  five 
times  and  the  baseline  airfoil  is  an  NACA  64A410. 

Shock-free  airfoils  obtained  at  the  same  flight  conditions  by 
varying  the  exponent  of  Equation  (lc)  and  thus  changing  the 
density's  dependence  on  flow  speed.  The  vertical  scale  is  mag- 
nified five  times  and  the  baseline  airfoil  is  an  NACA  64A410. 

Sonic  surface  for  the  shock-free  rectangular  wing  obtained  by 
modifying  a wing  with  a parabolic  arc  airfoil  section,  and  the 
pressure  coefficients  on  the  original  and  modified  wing  as 
calculated  by  the  numerical  algorithm  of  Ref.  19.  The  thickness 
distribution  of  the  baseline  wing  is  parabolic. 

Changes  required  in  the  surface  slope  at  various  lateral  sta- 
tions to  provide  shock-free  flow  over  Che  rectangular  wing  of 
Figure  9. 

Sonic  surface  on  the  shock-free  swept  wing  used  chat  corresponds 
to  the  design  pressure  coefficients  shown  in  Figure  12. 


Figure  12.  Comparison  of  the  computed  pressure  coefficient  on  a swept  wing, 
with  an  HACA  64A410  center  section  profile  and  an  elliptic 
thickness  distribution,  with  the  pressure  coefficient  obtained 
by  computing  the  flow  past  the  modified  wing  using  the  same 
numerical  algorithm.  The  leading  edge  sweep  is  30  DEG  and  the 
trailing  edge  sweep  15  DEG. 
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